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We consider some consequences of the presence of a cosmological lepton asymmetry in the form 
of neutrinos. A relic neutrino degeneracy enhances the contribution of massive neutrinos to the 
present energy density of the Universe, and modifies the power spectrum of radiation and matter. 
Comparing with current observations of cosmic microwave background anisotropies and large scale 
structure, we derive some constraints on the neutrino degeneracy and on the spectral index in the 
case of a flat Universe with a cosmological constant. 
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I. INTRODUCTION 



P-(' 

It is generally assumed that our Universe contains an approximately equal amount of leptons and antilcptons. The 
lepton asymmetry would be of the same order as the baryon asymmetry, which is very small as required by Big Bang 
Nucleosynthesis (BBN) considerations. The existence of a large lepton asymmetry is restricted to be in the form of 
neutrinos from the requirement of universal electric neutrality, and the possibility of a large neutrino asymmetry is 

■ still open. From a particle physics point of view, a lepton asymmetry can be generated by an Affleck-Dine mechanism 
|E| without producing a large baryon asymmetry (see ref. [Q for a recent model), or even by active-sterile neutrino 
oscillations after the electroweak phase transition (but in this last case, it might not be of order unity) H . This lepton 
asymmetry can postpone symmetry restoration in non-supersymmetric or supersymmetric models 0] (note that this 
is also true for other charges |J). 

In this paper we study some cosmological implications of relic degenerate neutrinos^. We do not consider any 
specific model for generating such an asymmetry, and just assume that it was created well before neutrinos decouple 

■ from the rest of the plasma. An asymmetry of order one or larger can have crucial effects on the global evolution of 
| the Universe. Among other effects, it changes the decoupling temperature of neutrinos, the primordial production of 

light elements at BBN, the time of equality between radiation and matter, or the contribution of relic neutrinos to 
the present energy density of the Universe. The latter changes affect the evolution of perturbations in the Universe. 
Qh. We focus on the anisotropies of the Cosmic Microwave Background (CMB), and on the distribution of Large Scale 
Structure (LSS). We calculate the power spectrum of both quantities, in the case of massless degenerate neutrinos, 
and also for neutrinos with a mass of 0.07 eV, as suggested to explain the experimental evidence of atmospheric 
neutrino oscillations at Super-Kamiokande Q . The cosmological implications of neutrinos with such a small mass are 
known to be very small, but we will see that this conclusion is modified if a large neutrino degeneracy exists. We also 
include in our analysis the possibility that the dominant contribution to the present energy density in the Universe 
is due to a cosmological constant: Ha ~ 0.6 — 0.7, keeping the Universe flat (fio + = 1), as suggested by recent 
observations (see [Q and references therein). 

The effect of neutrino degeneracy on the LSS power spectrum was studied in ref. as a way of improving the 
agreement with observations of mixed dark matter models with eV neutrinos, in the case of high values of the Hubble 
parameter. Also, Adams & Sarkar || calculated the CMB anisotropies and the matter power spectrum, and compared 
them with observations in the fl\ = case for massless degenerate neutrinos. More recently, Kinney & Riotto (Tc[ | 
also calculated the CMB anisotropies for massless degenerate neutrinos in the f^A = 0.7 case. 

This paper is organized as follows. In section H we calculate the contribution of massive degenerate neutrinos 



to the present energy density of the Universe. In section III, we explain how to calculate the power spectra, with 
the help of the code cmbf ast In section |v|, we discuss the effect of the degeneracy on CMB anisotropies and 
the matter power spectrum, both for massless neutrinos and m v = 0.07 eV. Finally, in section |v|, we derive some 
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1 Here, by degeneracy, we mean that there exists a large neutrino-antineutrino asymmetry, or vice versa, and not a degeneracy 
in the mass sense. 
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constraints on the neutrino degeneracy from CMB and LSS data, in the particular case of a flat Universe with an 
arbitrary cosmological constant and for standard values of other cosmological parameters. 



II. ENERGY DENSITY OF MASSIVE DEGENERATE NEUTRINOS 

The energy density of one species of massive degenerate neutrinos and antineutrinos, described by the distribution 
functions f„ and / p , is (here and throughout the paper we use Ti = c = ks = 1 units) 

<fp . , , . ,. 1 



Pu + p-u = J E v {f v (p) + Mp)) = ^J q dp p 2 Vp 2 + ml(f v (p) + Mp)) , (l) 

valid at any moment. Here p is the magnitude of the 3-momentum and m v is the neutrino mass. 

When the early Universe was hot enough, the neutrinos were in equilibrium with the rest of the plasma via the 
weak interactions. In that case the distribution functions /„ and fp changed with the Universe expansion, keeping 
the form of a Fermi-Dirac distribution, 

Mp) = 7 — - — \ Mp) = 7 — - — x (2) 

P M ,1 I P , M 



Here /j, is the neutrino chemical potential, which is nonzero if a neutrino-antineutrino asymmetry has been previously 
produced. Later the neutrinos decoupled when they were still relativistic^], and from that moment the neutrino 
momenta just changed according to the cosmological redshift. If a is the expansion factor of the Universe, the 
neutrino momentum decreases keeping ap constant. At the same time the neutrino degeneracy parameter £ = n/T v 
is conserved, with a value equal to that at the moment of decoupling. Therefore one can still calculate the energy 
density of neutrinos now from eq. (|]) and eq. (^), replacing fx/T v by £ and p/T v by pj (y v To), where To — 2.726 K and 
y u is the present ratio of neutrino and photon temperatures, which is not unity because once decoupled the neutrinos 
did not share the entropy transfer to photons from the successive particle annihilations that occurred in the early 
Universe. In the standard case, the massless non-degenerate neutrinos decoupled just before the electron-positron 
pairs annihilated to photons, from which one gets the standard factor y v = (4/11) 1 / 3 . 

In the presence of a significant neutrino degeneracy £ the decoupling temperature T(£) is higher than in the 
standard case |T^;|l3|. The reaction rate T of the weak processes, that keep the neutrinos in equilibrium with the 
other species, is reduced because some of the initial or final neutrino states will be occupied. The authors of ref. |l3[ | 
used the Boltzmann equation to calculate T for the process Vd + ^d e + + e~ (here Vd denotes degenerate neutrinos), 
including the corresponding Fermi blocking factors. It was found that the neutrino decoupling temperature is Td ec (S,) ~ 
0.2£ 2 / 3 exp(£/3) MeV (for or v T ). Therefore if £ is large enough, the degenerate neutrinos decouple before the 
temperature of the Universe drops below the different mass thresholds, and are not heated by the particle- antiparticle 
annihilations. The ratio of neutrino and photon temperatures is thus reduced accordingly. 

The present contribution of these degenerate neutrinos to the energy density of the Universe can be parametrized 
as p v = 10 4 /i 2 f2„ eV cm -3 , where f2„ is the neutrino energy density in units of the critical density p c = 3H 2 Mp/8TT, 
Mp = 1.22 x 10 19 GeV is the Planck mass and H — lOO/i Km s _1 Mpc -1 is the Hubble parameter. The value of 
p u can be calculated as a function of the neutrino mass and the neutrino degeneracy £, or equivalently the present 
neutrino asymmetry L v defined as the following ratio of number densities 

n 7 12C(3) 

We show^ in figures [l] and |^ the contours in the (m„ , L v ) and (m u , £) planes that correspond to some particular values 
of h 2 flp. One can see from the figures that there are two limits: massive non-degenerate neutrinos and massless 
degenerate neutrinos. The first case corresponds to the vertical lines, when one recovers the well-known bound on 
the neutrino mass m v < 46 eV for /i 2 f2„ = 0.5. On the other hand, for very light neutrinos, the horizontal lines 
set a maximum value on the neutrino degeneracy, that would correspond to a present neutrino chemical potential 



2 Unless the neutrino mass is comparable to the decoupling temperature, 0(m„) ~ MeV. 

3 Here we assume £ > 0, but the results are also valid for £ < provided that £ and L„ are understood as moduli. 
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Mo 5s 7.4 x 10 3 eV, also for h 2 £l v — 0.5. In the intermediate region of the figures the neutrino energy density is 
p v ~ Tii v n v (£) and the contours follow roughly the relations 

L v (^) ~ 24.2fc a n„ 



A similar calculation has been recently performed in reference |14 . However the difference between neutrino and 
photon temperatures was not properly taken into account for large £. It was argued that, since the number density 
of highly degenerate neutrinos is larger than in the non-degenerate case, the neutrinos would have been longer in 
thermal contact with e + e~, therefore sharing with photons the entropy release. However this is not the case |L3| as 
we discussed before. 

The presence of a neutrino degeneracy can modify the outcome of Big Bang Nucleosynthesis (for a review see ftsf ). 
First a larger neutrino energy density increases the expansion rate of the Universe, thus enhancing the primordial 
abundance of 4 He. This is valid for a nonzero £ of any neutrino flavor. In addition if the degenerate neutrinos are of 
electron type, they have a direct influence over the weak processes that interconvert neutrons and protons. This last 
effect depends on the sign of £„ e . Both effects may be simultaneously important and it could be possible in principle 
to explain the observed primordial abundances with a large baryon density, Vlsh? ~ 1 p^JI^ |. However this possibility 
is ruled out by the fact that in that case our Universe would have been radiation dominated during a longer period 
and the observed large-scale structure would be difficult to explain. From BBN one gets the following constraint 

- 0.06 < f„ e < 1.1 (5) 

while a sufficiently long matter dominated epoch requires 

l&wvl<6.9 (6) 

This estimate from fl3}| agrees with our analysis in section Assuming that the degenerate neutrinos are or v T , 
this places a limit on the degeneracy as shown by the horizontal line in figures [l] and §. 

III. POWER SPECTRA CALCULATION 

We compute the power spectra of CMB anisotropies and large-scale structure using the Boltzmann code cmbf ast 
by Seljak & Zaldarriaga |11| , adapted to the case of one family of degenerate neutrinos (v, D), with mass m v and 
degeneracy parameter £. Let us first review the required modifications. We use the notations of Ma & Bertschinger 
|L6|], and for all issues not specific to our case, we refer the reader to this review. 

Background quantities can be rewritten in terms of two dimensionless parameters (M, Q) 

_ mv_ m„(eV) = ap_ 

T uQ 8.6170 x 10- 5 x (4/ll)V3 T (K)' T M K > 

(we are assuming y v — (4/11) 1 / 3 , and therefore £ < 12 Jl3|]; the scale factor is defined so that a = 1 today). For Super- 
Kamiokande neutrinos with m„ = 0.07 eV, M ~ 417. We then get for the mean density, pressure and phase-space 
distributions 

1 V 

2tt 



P»+P» = 7T^ I Q 2 dQ\/Q 2 + anP (MQ) + UQ)) , 



n q 2 



6?r J VQ + a M 
MQ) = o-l : i ■ MQ) = ' 



In the case of massive degenerate neutrinos, these integrals must be calculated for each value of the scale factor, and 
also at the beginning of the code in order to find f2„ today. On the other hand, for massless neutrinos, there is an 
exact analytic solution 
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(9) 



So, if we define an effective number of massless neutrino families N e ff = 3 + 30/7(£/7r) 2 + 15/7(£/-7r) 4 , the mean 
density and pressure for all neutrinos will be given by these ones for one massless non-degenerate family, multiplied 
by N eff . 

Let us now consider perturbed quantities. We define *&„ and the perturbations of the phase space distribution 
for v and v, through 



Sf v (x, Q,n,r) = f v {Q)^ v {x,Q,n,T) 
6fu(x,Q,h,T) = fp(Q)^ i? (x,Q,n,T) 



(10) 



(n is the momentum direction: p = pn). For our purpose, which is to integrate the linearized Einstein equations, it 
can be shown that only the following linear combination is relevant 



fv^v + fu^u 

fu + fp 



(11) 



Using the Boltzmann equations for and \&p, it is straightforward to show that the evolution of '5 (in Fourier space 
and in the synchronous gauge, see [fill, eq. (40)) obeys 
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d\n{f v + U) 
d\nQ 



h + 6i) 



{k.nf 







(12) 



This equation depends on £ only through the last term, which is the gravitational source term. 

In the case £ = 0, the quantity (dln(/ l/ )/dlnQ) has a simple interpretation: it is the Q-dependence of a planckian 
perturbation of the phase space distribution. In other words, a shift of the blackbody temperature AT/T(i, n,r) 
corresponds to a perturbation 



(13) 



Since the gravitational source term in the Boltzmann equation is proportional to this quantity, the planckian shape 
is unaltered for massless neutrinos, and also for massive neutrinos when they are still relativistic (indeed, when 
Q 2 3> a 2 M 2 , the Q-dependence of the Boltzmann equation ( |l2|) vanishes in the second term, and remains only in the 
third term). When £ ^ 0, the source term in eq. (O) is proportional to 



rfln(/„ + / p ) 
dlnQ 



Q(l + ch£ chQ) 
(che + e-Q)(cK + chQ)' 



(14) 



When neutrinos are still relativistic, W is proportional to this quantity, even if it cannot be simply interpreted in 
terms of blackbody temperature perturbations. 

We can now specify all the changes required in cmbf ast, first in the case of massive degenerate neutrinos. As usual, 
^ can be expanded in a Legendre series: \& = Y^o(~ *) + l^i-Pf- It is easy to show that for each multipole 
\E';, the evolution equation and the initial condition are both identical to those of the non-degenerate case, provided 
that we replace (cZln(/ J/ )/dln(5) by eq. (|l4|). So, in summary, one only needs to modify the homogeneous phase- 
space distribution, its logarithmic derivative with respect to Q, and the initial calculation of f£„. Also, in order to 
obtain a good precision in the CMB anisotropy spectra, one must set I = 5 for the number of multipoles 'J; to be 
time-integrated. For transfer functions, the value I = 25 proposed by the code is sufficient. 

In the case of massless degenerate neutrinos, the situation is even simpler. The Q-dependence of the Boltzmann 
equation can be integrated away, just like in the non-degenerate case. For this purpose, we must introduce the 
Q-independcnt variable F v 



(15) 



1=0 



and integrate by part the last term in equation (|l2|). The multipoles F v i are exactly identical for degenerate and 
non-degenerate massless neutrinos, because they share the same evolution equations and initial conditions. So, the 
effect of £ arises only through the background quantities in eq. (|^) and is completely described by introducing an 
effective number of massless neutrinos. 
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IV. RESULTS 



First, as a consistency check, we compute CMB anisotropics and transfer functions for different values of £, choosing 
a very small mass m v < 0.001 eV. We check that the results match exactly those obtained with the unmodified version 
of cmbf ast, when the appropriate effective neutrino number N e ff is specified. 

The effect of £ and m„ on the CMB anisotropy spectrum can be seen in figure ||. We choose a set of cosmological 
parameters (h = 0.65, fib = 0.05, f^A = 0.70, flcDM = \ — fib — flu — fl\, Q rm s-ps = 18 /xK, flat primordial 
spectrum, no reionization, no tensor contribution), and we vary £ from to 5, both in the case of massless degenerate 
neutrinos (solid lines) and degenerate neutrinos with ?7i„ = 0.07 eV (dashed lines). Let us first comment the massless 
case. The main effect of £ is to boost the amplitude of the first peak^ Indeed, increasing the energy density of 
radiation delays matter-radiation equality, which is known to boost the acoustic peaks fl7j] (the same explanation 
holds for the effect of fl\ in flat models). For the same reason, all peaks are shifted to higher multipoles, by a factor 
((1 + a eq / a*) 1 / 2 — (agq/a*) 1 / 2 ) -1 ji~7| (a eq is the scale factor at equality, and increases with £, while the recombination 
scale factor a* is almost independent of the radiation energy density). Secondary peaks are then more affected by 
diffusion damping at large /, and their amplitude can decrease with £. 

In the case of degenerate neutrinos with m u — 0.07 eV, the results are quite similar in first approximation. Indeed, 
the effects described previously depend on the energy density of neutrinos at equality. At that time, they are still 
relativistic, and identical to massless neutrinos with equal degeneracy parameter. However, with a large degeneracy, 
flu today becomes significant: for £ = 5, one has fl v — 0.028, i.e. the same order of magnitude as fib- Since we 
are studying flat models, fl v — 0.028 must be compensated by less baryons, cold dark matter (CDM) or In our 
example, and fl\ are fixed, while flcDM slightly decreases. This explains the small enhancement of the first peak 
compared to the massless case. Even if this effect is indirect, it is nevertheless detectable in principle (even if one does 
not impose the flatness condition, the effect of fl v will be visible through a modification of the curvature). In figure ||, 
for £ = 0, the first peak maximum is enhanced by only 0.37%, while for £ = 5, there is an increase of 3.4%, detectable 
by the future satellite missions MAP and Planck, unless there arc large parameter degeneracies. It is well-known that 
such degeneracies are generally removed when CMB and LSS data are combined for parameter extraction Jl8| . 

We plot in figure | the power spectrum P(k) for the same models as in figure |3| normalized on large scales to 
COBE. The effect of both parameters £ and m„ is now to suppress the power on small scales. Indeed, increasing £ 
postpones matter-radiation equality, allowing less growth for fluctuations crossing the Hubble radius during radiation 
domination. Adding a small mass affects the recent evolution of fluctuations, and has now a direct effect: when the 
degenerate neutrinos become non-relativistic, their free-streaming suppresses growth of fluctuations for scales within 
the Hubble radius. For non-degenerate neutrinos, this effect is known to reduce power on those scales by a relative 
amount AP/P ~ 8fl„/flo [[l9| (we introduced flo = 1 — ^a)- So, even with m u = 0.07 eV and £ = 0, it is significant, 
especially at low fl . In the models of figure^, P(k) decreases by ~ 5 %, in agreement with the theoretical prediction 
(fl v = 1.8 x 10~ 3 , fio = 0.3). However, at £ = 5 (i.e. fl v — 0.028), this effect is even larger: P(k) decreases by a factor 
2.2, instead of an expected 1.7. This effect is likely related to the phase-space distribution of neutrinos with a chemical 
potential: their average momentum is shifted to larger values, making the free-streaming suppression mechanism even 
more efficient. 

Let us compare our results with those of previous works. The effect of £ on the CMB for massless neutrinos and 
fl\ = is the same as that one found in j^] . We also agree with the revised results in jl(| . 



V. COMPARISON WITH OBSERVATIONS 



Since the degeneracy increases dramatically the amplitude of the first CMB peak, we expect large £ values to be 
disfavored in the case of cosmological models known to predict a fairly high peak. On the other hand, a high £ is 
likely to be allowed (or even favored) for models that predict systematically a low peak, unless a large scalar spectral 
index n > 1.2 {blue tilt) is invoked. For instance, the degeneracy is likely to be favored by: (i) a large contribution 
of tensor perturbations; (ii) a significant effect from reionization; (iii) a low baryon density; (iv) a large h (h > 0.7); 
(v) flat models with fl\ < 0.6; etc. For such models, the peak amplitude can be boosted by £, keeping n close to 
one, which is more natural from the point of view of inflation. However, a careful case-by-case analysis is required, 



4 In fact, this is not true for very large values of £. In such cases, recombination can take place still at the end of radiation 
domination, and anisotropies are suppressed. For our choice of cosmological parameters, this happens for £ > 7, but in such a 
case the location of the first peak is I > 450, and the matter power spectrum is strongly suppressed. 
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since the effects of £ and n on CMB and LSS spectra are far from being equivalent. Our goal here is not to explore 
systematically all possibilities, but to briefly illustrate how £ can be constrained by current observations for flat models 
with different values of £!a- Recent results from supernovae [^0|, combined with CMB constraints, favor flat models 
with n A ~ 0.6- 0.7 §. 

We choose a flat model with h = 0.65, Qb = 0.05, Q r ms—ps = 18 /iK, no rcionization and no tensor contribution, 
and look for the allowed region in the space of free parameters (12a , £j n )- The allowed region will not be defined using 
a maximum likelihood analysis, but with the more conservative technique called "concordance" by Wang et al. pi) , 
which consists in taking the intersection of regions allowed by each experiment. 

For simplicity, we take into account only a few constraints on the matter power spectrum, known to be representative 
of the large amount of available data: the value of a% (the variance of mass fluctuations in a sphere of radius 
R = 8/i _1 Mpc) given for flat models in Q, at 95 % confidence level (CL)g; a \ 2 comparison with the STROMLO- 
APM redshift survey ^4j, at scales well within the linear regime, also with 95 % CL[| and finally, the constraint on 
bulk velocity at R = 50k" 1 Mpc p5[ , taking into account the cosmic variance. Except for the updated a$ constraint, 
we use exactly the same experimental tests as in 26 1, and refer the reader to this paper for details. For CMB data, 
we perform a % 2 analysis based on 19 experimental points and window functions, taking into account the Saskatoon 
calibration uncertainty, in the way suggested by [ p7| . The list of data that we use is given in pt| , and again allowed 
regions correspond to 95 % CLQ. We do not take into account the most recent experiments, for which window functions 
are still unpublished; they are anyway in good agreement with the data considered here. 

We plot in figure | the LSS and CMB allowed regions in (£, n) parameter space, corresponding to J7a = 0, 0.5, 0.6, 0.7. 
For Oa = 0.5 — 0.7, the LSS window just comes out of as limits. For Qa = 0, the lower LSS constraint is from as, 
and the upper one from APM data. In the case of degenerate neutrinos with m v = 0.07 eV, the windows are slightly 
shifted at large £, since, as we saw, the effect of £ is enhanced (dotted lines on the figure). The CMB allowed regions 
do not show this distinction, given the smallness of the effect and the imprecision of the data. One can immediately 
see that LSS and CMB constraints on n are shifted in opposite direction with £: indeed, the effects of £ and n both 
produce a higher CMB peak, while to a certain extent they compensate each other in P(k). So, for J7a = 0.7, a case 
in which a power spectrum normalized to both COBE and as yields a too high peak^, a neutrino degeneracy can only 
make things worst. On the other hand, for J7a = 0.5 — 0.6, a good agreement is found up to £ ~ 3. 

Let us finally consider the Q\ = case in which, after COBE normalization of the power spectrum, there is a 
well known discrepancy between the amplitude required by as and the shape probed by redshift surveys: these two 
constraints favor different values of n. We find that the neutrino degeneracy can solve this problem with £ > 3.5; 
however, the allowed window is cut at £ ~ 6 by CMB data, and we are left with an interesting region in which On = 1 
models are viable. This result is consistent with However, current evidences for a low Slo Universe f|(],|2£| are 
independent of the constraints used here, so there are not many motivations at the moment to consider this window 
seriously. 



VI. CONCLUSIONS 



We have considered some cosmological implications of a large relic neutrino degeneracy. We have shown that 
this degeneracy enhances the contribution of massive neutrinos to the present energy density of the Universe. For 
instance, neutrinos with a small mass m v ~ 10~ 2 eV can contribute significantly to f^Oi provided that there is a large 
neutrino-antineutrino asymmetry. 

Our main result is the computation of the power spectra of CMB anisotropies and matter density in presence of a 
neutrino degeneracy. We found, in agreement with H , that it boosts the amplitude of the first CMB peak, shifts the 
peaks to larger multipoles, and supresses small scale matter fluctuations. These effects follow the increase of neutrino 
energy density, that delays matter-radiation equality. 

We extended the calculation to the case of massive degenerate neutrinos, and showed the results for a mass of 0.07 
eV, as suggested by the Super-Kamiokande experiment. This mass has a small effect on CMB anisotropies. Indeed, 



5 The Viana & Liddle result JES is in very good agreement with an independent derivation by Girardi et al. Q. 

6 This confidence level stands for the goodness-of-fit of the model: when x 2 is greater than some value, the probability that 
we find the observed dataset, assuming the model to be valid, is smaller than 5 %. For the APM data, we have 9 — 3 degrees 
of freedom, and the limiting value is found in numerical tables to be x 2 — 12.5. 

7 Here we have 19-3 d.o.f.; the 95 % CL is given by \ 2 = 26. 

8 At least, for the values of the other cosmological parameters considered here. This situation can be easily improved, for 
instance, with h = 0.7. 
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such light neutrinos are still relativistic at recombination, but in presence of a degeneracy, they can account for a 
substantial part of the density today, of order f2„ ~ 10~ 2 . Also, we showed that small scale matter fluctuations are 
much more suppressed when the degenerate neutrinos are massive, because free-streaming of non-relativistic neutrinos 
is more efficient when their average momentum is boosted by the chemical potential. 

We compared our results with observations, in the restricted case of a flat universe with arbitrary (Oa,£, n) and 
fixed values of other cosmological parameters. We found that for J7a — 0.5 — 0.6, a large degeneracy is allowed, up 
to £ ~ 3. However, this upper bound is smaller than the value £ ~ 4.6 needed to explain the generation of ultra-high 
energy cosmic rays by the annihilation of high-energetic neutrinos on relic neutrinos with mass m„ = 0.07 eV [ |29| . 
We also tried smaller values of even if they are not favored by combined CMB and supernovae data. It turns out 
that a large degeneracy can account for both CMB and LSS constraints even with Qq = 1, provided that 3.5 < £ < 6. 
This analysis could be extended to other cosmological models. For instance, the degeneracy is likely to be compatible 
with a large contribution of tensor perturbations to large scale CMB anisotropies. 

Finally, it turns out that the degeneracy parameter and the mass of degenerate neutrinos have effects within the 
level of detectability of future CMB observations and redshift surveys, even with m v ~ 0.07 eV. However, a careful 
analysis should be performed in order to detect possible parameter degeneracy between £, m v and other cosmological 
parameters. 
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FIG. 1. Present energy density of massive degenerate neutrinos as a function of the degeneracy f. The curves correspond to 
different values of h 2 Q, v and the horizontal line is the upper bound from eq. (0). 



FIG. 2. Same as figure y for the neutrino asymmetry L v . 

FIG. 3. CMB anisotropy spectrum for different models with one family of degenerate neutrinos. Solid lines account for the 
case of massless degenerate neutrinos, and correspond, from bottom to top, to £ = 0, 3, 5. Dashed lines correspond to degenerate 
neutrinos with mass m v = 0.07 eV. Other parameters are fixed to h = 0.65, Q.b = 0.05, Qa = 0.70, flcDM — 1 — Hf, — — ^a, 
Qrms-ps = 18 nK, n = 1. We neglect reionization and tensor contribution. 

FIG. 4. Present power spectrum of matter density, for the same parameters as in the previous figure. From top to bottom, 
to £ = 0,3,5. 

FIG. 5. LSS and CMB constraints in (£, n) space, for different choices of Q.a'- from top left to bottom right, 
Oa = 0, 0.5, 0.6, 0.7. The underlying cosmological model is flat, with h = 0.65, Qb ~ 0.05, Q r ms- P s = 18 fjK, no reion- 
ization, no tensor contribution. The allowed regions are those where the labels are. For LSS constraints, we can distinguish 
between degenerate neutrinos with m v = (solid lines) and m„ = 0.07 eV (dotted lines). 
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